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ABSTRACT 


This thesis examines the problem of predicting distortion 
induced stall’ in an axial flow compressor. Three instantane- 
ous distortion indices are used on the data presented. The 
data consist of eleven stall events run at 2.5 Mach by NASA 
‘Lewis’-on a J85-GE-13 turbojet engine. The three indices 
used include two indices based on vorticity and developed at 
the Naval Postgraduate School. The third index was developed 
by Pratt and Whitney and is based on total pressure distri- 
butions. The three indices are compared by visual means us- 
ing computer plots and are also compared by using a simple 
statistical approach. The computer program used was first 
developed by Shoemaker in his Master's thesis. The program 
is expanded in this thesis to compute and compare the three 


indices. 
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I. INTRODUCTION 


This paper addresses the problem of distortion-induced 
Semoriit Stall in an axial flow compressor. The problem is 
a current one in industry and has not been wholly resolved. 
The present method of dealing with the problem of distortion 
is to operate axial flow compressors with a sufficient stall 
margin to avoid the problem altogether, but this is an expen- 
Sive procedure. Best engine efficiencies are obtained by 
operating close to the stall or surge line. By providing a 
large stall margin, we are in effect sacrificing performance 
via reduced efficiencies. The problem is especially germane 
femmilitary aircraft. The normal operating conditions of 
military aircraft include hard maneuvering, high "g'" flight, 
steam ingestion from aircraft carrier launches, and exhaust 
gas ingestion from forward firing weapons. These operating 
regimes cause a pressure, temperature, and gaseous molecular 
weight distortion of the flow entering the compressor section 
of the engine. 

Mmocter to deal with distortion, it is first necessary 
to find or define an adequate measure of it. A distortion 
index based on steady-state or long-time averaged values of 
distortion has been attempted. This method has met with some 
degree of success, but it is apparent that peak instantane- 
ous distortions of sufficient magnitude can induce stall when 
the time averaged distortion factor remains relatively low. 


tiememrrent approach and the approach of this paper is to 
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examine distortion as an unsteady phenomenon and to measure 
distortion over time intervals comparable with the character- 
istic times of an axial compressor. 

A computer program was first developed by Ens. Shoemaker 
[ref 3] utilizing an existing interpolation scheme developed 
at the Lewis Research Center of NASA. The program calculates 
the total pressure distribution just prior to the compressor 
section of a J&85-GE-13 turbojet engine as a function of time. 
The program also utilizes the theory of Fuhs, Farmer, and 
Myerson jret £, ref 2] to compute the circumferential and 
radial components of vorticity across the face of the com- 
pressor section as a function of time. The computer program 
was extended for this paper to calculate and statistically 
compare two existing distortion indices plus an additional 
index developed by Downs and studied in this thesis [ref 4]. 
This procedure was applied to eleven stall events recorded 


and analyzed by the NASA Lewis Research Center [ref 9]. 


steal, 





II. NATURE OF THE PROBLEM 


A. COMPRESSOR SURGE 

The line:-marked ''surge" in Pare uine 1 denotes the locus of 
unstable operation points for an axial compressor. The 
cause of compressor surge may be visualized in the following 
manner. At any given reference RPM (N//6) and blade radius, 
a rotor speed U is specified. Consider the velocity diagram 
in figure 2. A reduction in referred mass flow rate will 
ese a reduction in the axial component of velocity. Since 
the rotor speed is fixed, an increase in the relative flow 
angle 8 will occur. Present compressor design incorporates 
a fixed rotor geometry and thus the relative flow angle may 
only increase to a maximum value which is dependent on biade 
design. 

An insight into the unstable nature of the surge line 
may be made by comparing the characteristics of the compres- 
[wieeto the characteristics of a throttle. A throttle is set 
to produce a certain back or receiver pressure, given an in- 
let pressure. We may visualize a system consisting of a 
compressor producing a high inlet pressure followed by a 
throttle which dumps the high compressor air to ambient con- 
meerons (fig 3). A plot of the characteristics of the hypo- 
thetical throttle and compressor are included in figure 4. 
The compressor is represented at a fixed reference RPM. The 
dashed lines represent fixed throttle settings. Selecting a 


particular throttle setting determines two possible 
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eempressor points. Consider as a throttle setting line B, 
wun possible operating points P and Q. First consider point 
P, If the mass flow rate is perturbed to a higher mass flow 
rate, the throttle setting will demand a higher pressure 
ratio, whereas the compressor produces a lower pressure ra- 
flor in this case the combination of compressor and throttle 
will return to point P. Conversely a lower flow rate will 
cause the throttle to demand a lower pressure ratio whereas 
the compressor will produce a higher ratio. Again the sys- 
tem will be forced to return 6 porn » ye Consider now. porn 
Q@. A slight perturbation mass flow will produce problems. 
The system is free in this case to follow any perturbations 
and the system is unstable. Neutral stability is represent- 
ed by point R where the slope of the throttle curve matches 


the slope of the compressor curve. 


B. CLEAN STALL VERSUS DISTORTION INDUCED STALL 

The surge line in figure 1 denotes the conditions at 
mumen stall will occur if the inlet velocity profile is uni- 
Zorn proaqucing a uniform total pressure distribution. In 
feimtality, the velocity profile is not uniform and pressure 
perturbations may produce a local condition over a small re- 
gion of the compressor face which will cause stall over that 
region. This stalled region in turn produces a condition 
mievineas rotating Stall. In figure 5, if the conditions are 
such that blade 2 becomes stalled, the effective flow area 


between blades 2 and 3 is reduced. This reduction in flow 


oe 
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oe 


Rotating Stall 


FIGURE 5 


erea causes a flow diversion which tends to overload and 


Stall blade 3. Continuation of this process will produce 


mobatbing stall. it is noteworthy that rotating stall propa- 


gates in the opposite direction from rotor velocity at about 


half rotor speed. Rotating 
fully developed stall. 

It is apparent that the 
tions across the compressor 
Steady-state conditions are 
locity distribution in turn 


bution which is the form of 


stall may develop into surge or 


instantaneous velocity distribu- 
face may cause stall when the 
below the surge line. This ve- 
produces a total pressure distri- 


the data examined in this thesis. 
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Pei DESTORTION INDICES 


PweeevORITICITY VIA CROCCO'S THEOREM 

Non-dimensional vorticity components in both the circum- 
ferential and radial directions of a cylindrical coordinate 
system may be calculated using the theory of Fuhs, Farmer and 
ibuermsconebref 1,.ref.2]. <dstarting with Crocco's Theorem and 
uSing an order of magnitude analysis, it was possible to re- 


duce the expression for these parameters to: 


weet ok UM Or, Ct a SE eee 
8 eee, Os ao r year Us 00 2r' 30 
where: 
Og WD ,  U _ ft 
Pr 2a a Tip 


For a detailed explanation, consult reference 1. The calcu- 
lations are examined more briefly in reference 2. 

Analytical studies and data analysis have been performed 
in Ph attempt to relate the effects of total temperature and 
total pressure distortions on compressor performance [ref 7}. 
The vorticity approach is well suited to this problem, as 
both total temperature and total pressure may be taken into 
Mmecount in the calculations of vorticity. It would be bene- 
ictal to examine data which consist of both temperature and 
pressure patterns across the compressor face. Additional 
work in this area would be helpful. The data examined in 
this thesis, however, did not present a Significant tempera- 


ture gradient and the expressions for vorticity could be 
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reduced to: 


oe SR _ ea aP' 
8 “aU Giese? ac yPir'Us 09 


The task of this thesis was to determine whether an index 
based on vorticity would provide a meaningful prediction of 
compressor stall. Two such indices based on vorticity are 
examined in this aeons The first index was proposed by 
Shoemaker [ref 3] and is based on circumferential vorticity. 
A second index proposed by Downs [ref 4] and based on both 
radial vorticity and total pressure is examined. In order 
to provide a means of comparison, a third index based on to- 


tal pressure and taken from Zonars [ref 8] was also examined. 


B. SHOEMAKER INDEX 

James E. Shoemaker approached the problem of finding a 
Wertacity based index by first printing out maps of pressure, 
Sercumirerential vorticity and radial vorticity. In order to 
accomplish this goal he utilized a set of subroutines devel- 
oped at Wright Patterson Air Force Base, and altered them to 
meant Out the maps at times near stall and at times far from 
ZeeeeecOo provide a comparison. A copy of his computer pro- 
Eamets included at the end of his Master's thesis [ref 3]. 
mone this procedure, Shoemaker found a vorticity pattern 
that frequently occurred within twenty-five milliseconds of 
stall, and was not common at other times. The pattern was 
characterized by a ring of large negative circumferential 
Portlet ty at the midspan of the compressor face, and a large 


positive ring at the tip. He developed an index based on 
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this pattern which will be referred to as the Shoemaker in- 
dex. He divided the compressor face into three concentric 
regions. Region one at r' less than 0.5 was ignored. Region 
II was established as that area between radii of 0.55 and 
0.70. Region III was that area between radii of 0.8 and 
0.95. The Shoemaker index then summed over Region II but set 
any positive circumferential vorticity equal to zero. Like- 
wise Region III was summed where any negative vorticity was 
set equal to zero. The instantaneous Shoemaker index was 


set equal to the sum of the absolute value of the two summa- 


meoms. in notation form, 
Shoemaker Index (T) = = Z tay Ne: > [w,] 
N 0 N 6 
a8 hee ; 
where Wy = positive values of circumferential vorticity in 
Region III 
W = negative values of circumferential vorticity in 
Region II 


In order to take into account the added effect of the 
duration of the pattern, Shoemaker set his final index equal 
to the product of three consecutive instantaneous indices de- 
scribed above. This latter procedure was not used in the 
present thesis. The Shoemaker index as used herein thus re- 
fers to the instantaneous index and is an indication of dis- 
Pert iOn in the radial direction. 

The Shoemaker index has one real advantage when the 
transport time from the measuring probes to the compressor 


Swaee (Where Stall originates) is considered. The pattern 
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develops at a time before stall which is larger than the 
transport time from the measuring probes to the stall region. 
mesoecould be very useful if variable stators or the fuel 


=omemoen were available to avert stall. 


C. DOWNS INDEX 

A second index which was proposed by Downs [ref 4] was 
formulated and is explored in this thesis. This index incor- 
porates both total pressure and radial vorticity. The basic 
idea was to consider the contribution of both total pressure 
and vorticity to the relative flow angle g. Increasing 8 
would cause the rotor blade incidence angle to increase and 
would tend to load the blade. A perturbation plot of §8 ver- 
sus 8, where the deviation of 8 from the reference line takes 
into account pressure and vorticity effects, is shown in fig- 
ure 6. The corresponding plot of total pressure and Ws is 
also shown. The reference line represents uniform flow. It 
can be seen that although W 1s dependent upon the pressure 
gradient, the combined effect of pressure and W, on 8 is de- 
pendent on both the gradient and relative magnitude of pres- 
sure. The Downs index was formulated to search the compressor 


face for the largest value of the combined effect of pressure 


plus W,, Over any 60 degree sector. In notation form, 
Bao 
Peveowtndex (t) = 2 [=——=—_—— + w') 
p max 60 
i 
where Pa = average total pressure over entire face of 


compressor at time t 
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Mires Wowns index 1S primarily an indication of distortion in 
the circumferential direction. For a more complete analysis, 


consult Downs [ref 4]. 


D. PRATT AND WHITNEY Kou DISTORTION INDEX 

In order to provide a means of comparison between vortic- 
meveebased indices and an index calculated directly from total 
pressures, the Pratt and Whitney Ko index was selected from 
Zonars [ref 8]. This index was selected primarily because 
of case four, which was one of the eleven cases investigated. 
Distortion indices are plotted for all eleven cases in Appen- 


mies, Neither of the vorticity indices had satisfactorily 


predicted stall in case four; from examining the pressure 


maps, Kaas appeared to be an appropriate index. The KOM index 
was applied to case four and to all other cases. This index 


was developed by Pratt and Whitney for the TF-30 installed in 


em@eer—-11l1A aircraft and is defined as follows: 


" See: Cian _ 
ky (eR), oc, 
ieee 
a ave et 
DM y om 


where c ratio of compressor inlet radius to ring radius 


1 = number of ring 
0 = largest continuous are of the ring over which 
the total pressure is below the ring average 
pressure 
je = ring maximum total pressure 
Tmax 
is = ring average total pressure 
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Pp : — Tine maninum total pressure 
Tmin 5 yy 


Mmeeean be Seen that the Kou index iS an indication of dis- 
tortion in the circumferential direction and puts more empha- 
Sis on the region near the hub. It would be expected that 
the Ko index and the Downs index would be comparable since 
wemesearch for circumferential distortion. This was in 


mee tchie Case and will be discussed in a later section. 


ZS 





IV. EXPLANATION OF DATA 


The NASA Lewis Research Center conducted an experimental 
wind-tunnel investigation to determine the effects of time- 
varying total-pressure distortions produced in a supersonic 
mete On a J&60-GE-13 turbojet engine {ref 5]. Selected data 
taken during this investigation were digitized on magnetic 
eae the) latter is the source of data used in this thesis. 
The following explanation of data is largely a brief summary 
of the explanation presented by Burstadt and Calogeras [ref 
5] and is given here to provide an overall understanding of 
the data analysis. 

The inlet was an axisymmetric mixed-compression type. 
The compressor was an eight-stage axial flow compressor with 
interstage bleed and variable inlet guide vanes. An after- 
burner section was present but was not used during the test 
procedure. The engine bleed system was modified to make the 
compressor operate with less than normal stall margin. The 
following excerpt is quoted from reference 5B: 

The test procedure was to establish an engine 

Seonarine® COndition in the vicinity of the com- 

pressor stall point. The FM tape recorder was 

started and a steady-state data scan was made. 

Then, while the FM tape recorder was still run- 

ning, the compressor pressure ratio was increased 

Slightly and the throttle adjusted to maintain a 

Seomoatant rotor speed. If stall did not occur, 

this procedure was repeated until the stall event 

was recorded on the FM tape. In Some cases the 

Stall occurred after the engine had stabilized 


at an equilibrium operating condition. These 
Webewecalled ‘drift' stalls. 





All eleven of the stalls examined in this thesis were of 
the "drift" type. NASA Lewis considered only seven of the 
eleven cases in their analysis. Cases 6, 8, 9, and 10 were 
excluded on the grounds either that their causes were attri- 
butable to steady-state distortion levels or that they showed 
very little loss in stall margin. Table I shows steady-state 
operating conditions for all eleven cases presented in this 
thesis. The NASA Lewis analysis identifies each case number 
as a steady-state reading number. These NASA identifying 
numbers are included in Table I for comparison [ref 5]. 

Steady-state and dynamic pressure instrumentation at the 
Sempressor face is shown in figure 7. Only the fluctuating 
component of pressure was recorded from each of the dynamic 
probes. Absolute pressures were later obtained by adding 
the fluctuating component to the steady-state value. 

Frequency response of the dynamic probes was flat to 
meoutr ZO00 hertz. The fluctuating signal from each pressure 
transducer was recorded on FM multiplexed tape at 152.4 cen- 
timeters per second. The data in a 225 millisecond period 
wast prior to stall were digitized at an effective rate of 
8000 points per second per channel. These data when processed 
by the computer program listed in the back of this thesis 
provided the distribution of total presSure across the com- 
PieescsOr face with a time spacing of 0.125 milliseconds between 
distributions. The distributions of total pressure were fur- 
Pact eprocessea to provide distributions of both circumferen- 


PE anCmec ati vorticity as a function of time. In order to 
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meauce the amount of calculating involved, only one out of 
je eereadings was actually processed so that there was five 
times 0.125 or 0.625 milliseconds of time between readings 
actually processed. 

The determination of stall time was made by NASA Lewis. 
The signals from the compressor Static pressure transducers 
were put on strip charts and examined by NASA. oe 
cases, the row of instrumentation at the 290 degree location 
emowed the earliest indication of stall. Since the origin 
of stall in a compressor stage causes an abrupt loss in the 


airflow pumping of that stage, a compression wave is 
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propagated upstream. This compression wave is referred to 
as a hammershock. The following excerpt is quoted from ref- 
erence 5D: 


Because there were only two axial rows of in- 
strumentation inside the compressor, it was not 
possible to determine exactly when stall originated 
inside the compressor. Although the first indica- 
tion was seen at the 290 degree location, a stall 
Zone could originate just past the 45 degree loca- 
tion but not be detected until it rotated past the 
instrumentation at 290 degrees. Since the stall 
zone rotated at about half rotor speed (8 msec/rev) 
there could be as much as a 6-millisecond delay be- 
fore it was detected at the 290 degree location. 


A one-millisecond period for transport time between sensors 
and stall stage, plus time for blade response, was predicted 
by NASA. This produces a six-pluS-one or seven-millisecond 
possible time lag between actual initialization of stall and 


fmm@erstall time listed in this thesis. 
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Vo wewunkibe INDICES APPLIED TO DATA 


A. EXAMINATION PROCEDURE 

The computer program first developed by Shoemaker [ref 
3) was adapted to compute each index as a function of time 
Hormeeche entire time interval provided by NASA Lewis. This 
time .interval vanied.slighily but was generally about 215 
Pameersecconds of data prior to first indication of stall. 

The DRAW subroutine of the Naval Postgraduate School was used 
to plot each index. These plots are included in Appendix A 
and are plots of the instantaneous version of each index. A 
numerical printout was also provided by the program. A sam- 
ple numerical printout for the Shoemaker index case one is 
included at the end of this thesis. Table II provides an 
Seeanation of the terms for the numerical printout. 

In addition to the computer plots, which may be used to 
determine if the given index is an accurate predictor of 
fwerl, a simple statistical approach was used. After the 
first fifty values of the index were computed, the mean and 
Standard deviation for those fifty values were obtained. 
Then, as the next fifty values of the index were computed, 
the following test was made. If the average of the three 
previous values of the index at any time eee i Certain 
@riterion value, the computer program printed out a stall 
warning with the time in milliseconds of occurrence. The 
criterion level was set in the following manner: 


ese emean + 2.505% Standard—deviation. 
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TABLE II 


EXPLANATION OF PROGRAM TERMS 


N = index number (a function of time). 

maeai(N) = value of instantaneous index at time "N". 

DFACT(N) = average of three previous FACT(N) values. 

TIME JetNet isecOonds. “OoLartsS at zero Lor 
Sec ele Wc a see 

PTAVE = average total pressure over compressor face at 
tames NN . 

MEAN = mean value of FACT(N) for 50 previous values of 
N, when not zero. 

STD-DEV = standard deviation of FACT(N) for 50 previous 


values of N, when not zero. 


At the end of each 50 values of the index, a new mean and 
standard-deviation were computed and the test criterion was 
established for the next 50 values of the index. This simple 
test was used on all thirty-three runs made, and the results 
are Summarized in the next section. Obviously, the test 
criterion could be lowered, with resulting greater occurrence 
of stall prediction, but this would also increase the number 
of false warnings. The value of test was initially set to 


Pecaret Stall for the K index Case three without any false 


DM 
wmairnanes. The test criterion could have been set to an opti- 
mumevalwe for each particular index with possible better re- 


Pwr cebitethig method was not used. The same test was applied 


pO wt | cGASsesS. Tiere tbeulcate method has the real advantage 
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meno being subject to interpretation. All thirty-three 
Sompucer plots are included in this thesis, since they cer- 
wainiy are subject to personal interpretation as to the pre- 


frerion Of compressor stall. 


B. COMPARISON OF RESULTS 

miere are two indicators of the results of the data re- 
M@meetron. First are the computer plots, which provide an 
eyeball examination of results. <A matrix summation of the 
computer plot results is provided by Table III. The second 
method was the use of statistical mean and standard-deviation 
discussed in the preceding section. A matrix summation of 
the statistical method is provided by Table IV. In both 
tables, where stall was predicted, the time to stall is given 
in milliseconds. 

There are several points of interest. As had been pre- 
viously expected, there is some correlation between the BoM 
index and the Downs iMidex. since HOUmerely on Circumierential 
fmPisetortion. Using the statistical method, all of the cases 


predicted by the K index were also predicted by the Downs 


DM 
index. The Downs index predicted two additional cases uti- 
ezine this method. Using the computer plot interpretation, 
the Kans method predicts six cases and the Downs index pre- 
mmers five, with four predictions common to both indices. 
Mieroiece were to use a combination of both vorticity based 


indices, the statistical method would have predicted eight 


Out of the eleven cases presented. Using this method there 


ll 
















Poise Lit 
COMPUTER PLOT SUMMARY 


(based on visual interpretation) 





RUN PRATT & WHITYEY  SHOEMAKSR DOWNS 
Ko, INDEX INDEX INDEX 









No Prediction Pe.5 Millisecs PFO nn Lilvcees 
















apo mrllisecs 25.0 millisecs No Prediction 
















4 same magnitude 


5 peo maliiscecs 


one false peak 
yO ili aceces 





Pl.Owmptlisecs 






4 false peaks 2 false peaks 








































4. Zeon ls secs No Prediction No Prediction 
a one false peak 
» Wom scaicvion 3,0 millisceos 
Oo smaller peaks 
6 No Prediction No Prediction Prediction 
T- No Prediction No Prediction Prediction 
8 millisecs No Prediction POs rysecs 





16.0 millisecs 2.0 millisecs 






millisecs 





false peak 5 same magnitude 





4.O millisecs No Prediction 2.0 millisecs 


No Prediction 


ee a a ee oe ome 
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| 6.0 millisecs No Prediction 


Oo larger peaxs 


totals | 6 Predictions 6 Predictions S Predictions 
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TABLE IV 
STATISTICAL SUMMARY 


RUN PRATT & WHITNEY SHOEMAKER DOWNS 


eam e INDEX INDEX 


No Prediction eo miliisecs Pen ecwnin Licsecs 

















1 false warning 2 false warnings 


No Prediction 25.9 millisecs No Prediction 


2 false warnings 


O.9 millisecs Soom iaSsecs 





BATE 





millisecs 


© false warnings 
NO. Predict2en 


false warning 


No Prediction No Prediction 


false warning 1 false warning 


Romtredi ction Ereovcotlon Noa Predaetaon 


fom ered 16.10) 


Prediction No Prediction 
7 : 





Prediction o Prediction 19.88 millisecs 















1 false warning 


2.88 millisecs Peeoo owt lisecs 





No Prediction 





2 false warninss 
f5,0 Mad bisecs 





Oe millisecs 9.0 milliseecs 





e false warnings 






1.62 millisecs No Prediction Peoaei dt secs 





{ 
No Prediction ious ecs No Prediction 


false warning 1 false warning 


totals ede pens WMO ePredictions OP ited vetlions 
Seen | 








would have been seventeen false warnings. The number of 
Pose warnings could be reduced by a reduction in the test 
@iemeersOn to a point where no loss in stall prediction would 
Secur. 

There could be as much as seven milliseconds of time lag 
between actual stall time and predicted stall time. This 
Meesents a problem since the data could include the effects 
of the hammershock which occurs after stall, if the time lag 
/mmeoo great. It is certainly no problem to produce an index 
which will sense the hammershock and the question arises, did 
maerwandex actually predict stall or were the effects of the 
hammershock being felt? By examining the numerical printout 
of average total pressure as a function of time and taking 
into account the possible seven-millisecond lag time, it was 
merarent that this could have been a problem in only four 
cases. These four cases should be held suspect and are cases 
3, $, 9, and 1i0. It is noteworthy that the Shoemaker index 
Peeeorcred Stall in two of these four cases sufficiently far 


Mmaeaavance to preclude the problem. 
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VI. CONCLUSIONS 


An investigation into the problem of predicting compres- 
sor stall using instantaneous distortion parameters was made. 
Some success was evidenced by both the statistical method 
paeecOomputer plot interpretation, but it was not possible to 
Peaick Stall in every case uSing one or any combination of 
the indices explored. This seems to indicate that an axial 
compressor is subject to many instantaneous distortion pat- 
terns as well as to steady-state distortion, whereas one in- 
dex may only search for one category of distortion. If the 
cause of compressor stall was 1n every case instantaneous 
meas GAistortion over the outer portion of the compressor 
face, then in every case the Shoemaker index should predict 
stall. If the cause was in every case circumferential dis- 
trotion over some arbitrary sixty degree sector then the 
Downs index would be successful. If the distortion conforms 
memerne ring pattern explored by the Pratt and Whitney index, 
then that method would succeed. In fact, it seems naive to 
believe that compressor stall could be predicted in every 
case uSing any one index. The indices explored in this the- 
sis and the indices used by industry today may certainly be 
improved upon, but it is doubtful that any one index will 
prove to be a panacea. 

In spite of the compound nature of compressor stall, in- 
stantaneous indices are useful. The Pratt and Whitney index 


used in this thesis was developed for the F-111A with a 





Seeciitiic inlet-engine combination. It is in that context 
that distortion parameters would appear most useful - that 
1S, to examine an inlet-engine system for compatibility. 
The performance of an inlet-engine system which too often 
experiences one type of distortion pattern may be improved 
Mom, GiStortion indices would be very useful in isolating a 
predominant pattern. 

The use of vorticity as an index parameter has some real 
advantages. Both temperature and pressure gradients may be 


macesrated into a single parameter. Vorticity may be visual- 


maecemore Simply and is a natural characteristic of the fluid. 


The Shoemaker index, when successful, harbingers stall suffi- 
S@entiy far in advance possibly to enable correction of the 

Gere rency UuSing engine controls. Vorticity appears to he a 
promising dimensionless parameter and deserves further inves- 


ie acion. 
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APPENDIX A 


COMPUTER. PLOTS 


The following thirty-three computer plots provide the 
instantaneous version of each of the three indices considered. 
In all cases, the horizontal scale is read in milliseconds. 
The vertical scale gives the index value and is a non-dimen- 
Sional number. No consideration was given to index magnitude 
a Dredicting stall. For any given case, a strong peak just 
Memere stall in relation to that particular case was consid- 


ered to be sufficient for stall prediction. 
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FIGURE A-2 
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FIGURE A-6 
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